Unusual Structure and Magnetism in MnO Nanoclusters 
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We report an unusual evolution of structure and magnetism in stoichiometric MnO clusters based 
on an extensive and unbiased search through the potential energy surface within density functional 
theory. The smaller clusters, containing up to five MnO units, adopt two-dimensional structures; 
and regardless of the size of the cluster, magnetic coupling is found to be antiferromagnetic in 
contrast to previous theoretical findings. Predicted structure and magnetism are strikingly different 
from the magnetic core of Mn-based molecular magnets, whereas they were previously argued to 
be similar. Both of these features are explained through the inherent electronic structures of the 
clusters. 



PACS numbers: 36.40. Cg, 75.75.-c, 71.15.Mb 

Transition metal oxide, especially MnO^^, clusters 
have recently attracted extensive multidisciplinary re- 
search activity because of their diverse and tunable mag- 
netic and catalytic properties. Generally, as compared 
to the bulk, the local magnetic moment is enhanced in 
smaller dimensions due to reduction in the number of 
neighboring atoms. This results in either an overall en- 
hancement of the total moment for the ferromagnetic 
(FM) case or lead to a finite moment even for the an- 
tiferromagnetic (AFM) case due to unequal compensa- 
tion of spin up and down electrons. Magnetic coupling 
also evolves with particle size, and such size evolution 
for MnO clusters is non-monotonic. MnO clusters with 
a diameter of 5-10 nm show FM behavioi"^ even though 
the bulk phase is AFM^. In contrast, Mn-based single 
molecular magnets (with magnetic core < 1.5 nm) show a 
'layered' AFM/ferrimagnetic structure within the mixed- 
valent Mn centers, resulting in a large magnetic moment 
and spin anisotrop}ffl^. Moreover, the MnO clusters take 
essential part in a variety of biological (catalytic) pro- 
cesses from photosynthesis to bacterially mediated or- 
ganic matter decomposition. The active inorganic cen- 
ter of the oxygen evolving photosystem II contains a 
manganese oxide cluster (Mn404Ca), which catalyzes the 
light-driven oxidation of watei!^. Indeed, synthetic com- 
plexes containing cuboidal Mn404 cores have been found 
to exhibit unique reactivity in water oxidation/02 evo- 
lution SI. 

The prediction of geometry at the atomic level is one of 
the most fundamental challenges in condensed matter sci- 
ence. The magnetic and catalytic properties (i.e., broadly 
speaking: the electronic structure) are strongly coupled 
to the 'inherent structure' (corresponding to minima of 
the potential energy surface (PES)) of the cluster. Ex- 
perimental evidence of structure and magnetic coupling, 
and their size evolution for the transition metal oxide 



clusters in the gas phase are scarce. However, the struc- 
ture and (FM) magnetism of (MnO)^^ clusters have been 
predicted theoreticall}^^^. Such a theoretical prediction 
is complex, and requires a systematic and rigorous search 
through the PES. This is essential to predict the deep- 
est minima. The complexity of the PES search increases 
with increasing cluster size. Possible geometrical struc- 
tures increase exponentially with cluster size, and for a 
given geometrical structure containing N magnetic ions, 
there are 2^ spin configurations (which may or may not 
be reduced depending on the symmetry). In contrast, 
all the previous theoretical attempts have been largely 
biased^ because: (i) geometrical structures were re- 
stricted by a particular symmetry, and (ii) magnetic 
structures were restricted to the FM regime. Such limited 
considerations search only a small subspace of the entire 
PES and thus, previously reported geometric/magnetic 
structures may not represent the true 'ground states'. In- 
deed, in contrast to the previous theoretical reports, in 
this communication we shall report, based on a rigorous 
PES search, that the (MnO)ic clusters show AFM cou- 
pling and also show unusual two-dimensional (2D) struc- 
tures up to a certain size. 

The spin-polarized density functional theory calcu- 
lations were conducted using; the VASP cod^ with 
the Perdew-Burke-Ernzerhof exchange-correlation func- 
tional^^ and the projector augmented wave pseudopo- 
tential^^ at a energy cut-off of 270 eV. Simple cubic su- 
percells were used with periodic boundary conditions, 
and it was made sure that that two neighboring im- 
age clusters were separated by at least 10 A of vacuum 
space. This ensured that the interaction of the clus- 
ter with its periodic image was negligible and reciprocal 
space integrations were carried out at the F point. We 
started with high symmetry structures, for example, a cu- 
bic structure (core of the Mni2-molecular magnet^) for 
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FIG. 1. (color online) Geometrical evolution shows a 2D to 
3D transition at (MnO)6. Minimum energy structures are 
always found to be antiferromagnetic. The up (down) Mn- 
atoms are shown with yellow (green) color. Oxygen atoms 
are shown in red. The Mn-0 (1.80 - 2.18 A) and Mn-Mn 
(2.55 - 3.10 A) distances, and the Mn-O-Mn angles (78°-108°) 
in these clusters are comparable to the Mn-based molecular 
magnet^^. 
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FIG. 2. (color online) (a) The trend in binding energy shows 
that the clusters adopt 2D structures until they contain five 
MnO units, above which they are 3D. The Mn-Mn magnetic 
coupling is always antiferromagnetic for the most stable so- 
lution. Moreover, the AFM coupling is always favorable over 
the FM coupling, in both 2D and 3D. (b) Total hybridization 
index, H, = 1-isp + 1-Lpd + Hsd, shows a direct correspondence 
to the energy. 



the (MnO)4 cluster. Spin-polarized Born-Oppenheimer 
molecular dynamics (BOMD) simulations at 1200 K were 
done for 20-30 ps to search for the lowest energy isomers. 
This approach could efficiently explore the PES. Several 
minimum energy structures were then picked from these 
BOMD simulations and were carefully restudied. All of 
these minimum energy structures were further optimized 
{local ionic relaxation) with all possible spin multiplic- 
ities. Moreover, we have also considered different spin 
arrangements for different atoms in the cluster for a par- 
ticular spin multiplicity. 

Although the magnetic structure is strongly coupled 
with cluster geometry, we first discuss the geometric 
evolution alone. For the MnO dimer, as expected, the 
Mn-0 distance is much smaller (1.65 A) than bulk, in 
agreement with experimental and diffusion Quantum 
Monte Carlo (DMC) results^^ Calculated Mn-0 stretch- 
ing frequency (920 cm~^) is slightly higher than the 
experimentally obtained values (832 cm~^ Ref.^^ and 
899 cm~^ Ref.^^) in the gas phase. The average Mn-0 
distances increase with cluster size, and the lowest en- 
ergy structures are shown in Fig. [T] Interestingly, the 
presence of oxygen alters the geometry of MnO clus- 
ter as compared to pure Mn clusters^^. Moreover, we 
find that clusters containing up to five MnO units ex- 
hibit two-dimensional (2D) structures (Fig. [T]), and the 
lowest-lying 3D structures are separated from these by 
a large energy difference [Fig. [2|a)]. The structural en- 
ergy difference A£^2D-3D for (MnO)4 and (MnO) 5 clus- 
ters are 1.28 and 1.11 eV (equivalent to internal vibra- 
tional temperatures of 1650 and 1075 K), respectively. 
However, the scenario is reversed from the (MnO)6 clus- 
ter (Fig. [!}, for which A£;3d-2D is 0.61 eV (472 K). 
The present findings contradict previous theoretical re- 
ports, where the PES search were highly biased^ ^. In 
contrast, DMC calculations biased with only FM cou- 



pling predicts a similar geometrical trend for x < 
The present results are in accordance with the predic- 
tion of mass spectra^^, where Ziemann and Castleman 
proposed non-cubic structural growth, and larger clus- 
ters were composed of relatively more stable hexagonal 
(MnO) 3 and rhombic (MnO) 2 units. Similarly, we also 
find that the (MnO)3 cluster is more stable (magic clus- 
ter) and serves as the building block for larger clusters. 
It was believed earlier that the small MnO clusters serve 
as the magnetic core of Mn-based single molecular mag- 
nets the present study confirms that there is no such 
structural resemblance in terms of symmetry. The pla- 
nar structure of (MnO)4 (Fig.[T]) is substantially different 
from the magnetic core of [Mn403Cl4(02CEt)3(py)3]2^, 
[Mni20i2(CH3COO)i6(H20)4]^ molecular magnets and 
the oxygen evolving center (Mn404Ca) in photosystem II 
^ . The structural symmetry of (MnO)6 is very different 
from that observed in octahedral [Mn6(/i6-0)(/i-OR)i2] 
molecular comple:xE^. However, the structural parame- 
ters such as Mn-0 (1.80-2.18 A), Mn-Mn (2.55-3.10 A) 
bond lengths, and Mn-O-Mn angles (78°-108°) in these 
clusters are similar as compared to the molecular magnets 
(1.85-2.23 A, 2.77-3.44 A, and 94°-137°, respectively )'ll^. 
It should also be noted that the Mn-0 distance in these 
clusters already approaches that of the bulk MnO (2.25 

A). 

The binding increases substantially due to the presence 
of oxygen (Fig. [2|. This can be explained in terms of 
electronic configurations. The free Mn atom has 3(i^45^ 
electronic configuration, and thus the Mn-dimer is weakly 
bonded^ ^. In contrast, our Bader charge analysiJ^ shows 
that Mn atoms lose 1.20e charge to O in (MnO) 2, which 
makes binding in these clusters stronger. This amount 
(~1.20e) of charge transfer remains the rule of thumb for 
all the minimum energy structures for the stoichiometric 
MnO clusters. In addition to the covalency, this signifi- 
cant charge transfer suggests an ionic contribution to the 
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Mn-0 bonding. 

Calculated stability, as measured by AEx = 
£;(MnO)3,_i +£;(MnO)^+i -2£;(MnO)^, where E is the 
total binding energy, shows local peaks for clusters with 
X = 3, 4 and 6, in agreement with experiments^. This 
indicate exceptional stability referring to their magic na- 
ture. Clusters with three and four MnO units serve as 
the building blocks for larger clusters. Two (MnO)3 clus- 
ter units are stacked in 3D for the ground state (Fig. [T]), 
when the similar stacking in 2D (with similar AFM mag- 
netic ordering) is 0.61 eV higher in energy [Fig. [2|a)]. 
Similarly, as shown in Fig.jl] a (MnO)3 unit is stacked in 
3D with a distorted (MnO) 4 unit to form the most stable 
(MnO)7 cluster. Finally, the (MnO)8 cluster is formed by 
staking two (MnO)4 units. 

To understand this morphological transition, we plot 
the orbital projected density of states (DOS) for the clus- 
ters (Fig. [3|. These show a clear trend: the energy 
spread of the orbitals is higher for the minimum energy 
structure, and also has a larger orbital overlap. For ex- 
ample, a (MnO)4 cluster has larger orbital spread and 
also has higher s-d and p-d hybridization in 2D compared 
to the respective optimal 3D structure. In contrast, the 
situation is reversed for (MnO)6 cluster, where the larger 
spread and higher s-d and p-d hybridization makes the 
3D structure favorable over the 2D one. The orbital hy- 
bridization can be quantified, and could explain the clus- 
ter morphology.t^ We calculate the {k-l) hybridization 
index, 

nki = Y.T.'^ikyju, (1) 

/ i 

where /c, / are the orbital indices, w^j^{w(i) is the square 
projection of the i-th Kohn-Sham orbital on to the k (/) 
spherical harmonics centered at atom / and integrated 
over an atomic sphere (radius depends on the atom type, 
Mn/CP^. Note that the spin index is inherent. How- 
ever, unlike gold clusters, in a system with active 0-p 
electrons, ^s in addition to the s-d hybridization index 
{^sd)i the 1-Lpd and 1-Lsp would also play an important 
role for determining the dimensionality. Indeed, we find 
[Fig. [2|b)] that the total hybridization index 1-L is always 
higher for the lowest energy structures for the entire size 
range. 

Next we turn our attention to the magnetic ordering. 
We find that the Mn-atoms in all clusters are AFM cou- 
pled. This gives rise to a net or 5 jiB moment for the 
clusters with even and odd number of MnO units, respec- 
tively. It is also interesting to note here that regardless 
of their dimensional nature, the AFM arrangement is al- 
ways energetically favorable [Fig.[2|a)]. Experimental ev- 
idences of magnetic structure in such clusters are scarce, 
but the calculated results for the MnO dimer (sextet) is 
in agreement with the only available experiment Al- 
though the moment is localized on the Mn-sites, a small 
p polarization (0.02-0.15 /j^b) '^^ observed for the oxygen 
atoms. These results of the magnetic structure are in 
direct contradiction with the previous results including 




FIG. 3. (color online) Orbital projected density states 
summed over all the atoms in the cluster (for 2D and 3D 
(MnO)4 and (MnO)6 clusters in their AFM state; FM state 
for 3D (MnO) 6 is also shown) show that the minimum en- 
ergy structures (both geometric and magnetic) correspond to 
larger energy spread of the orbitals and higher hybridization. 
Gaussian smearing (0.1 eV) has been used. 



the DMC calculations^ As mentioned before, all the 
previous studies were highly biased and only scanned a 
small subspace of the PES constrained either by geomet- 
rical structure with high symmetry or by FM coupling, or 
constrained by both. We find that the presence of oxygen 
stabilizes Mn-Mn AFM coupling. For example, the en- 
ergy difference A£:(FM - AFM) is only -0.2 and 0.14 eV 
for pure Mn4 and Mne clusters respect ivel}^^, whereas 
this difference substantially increases to 0.92 and 0.67 
eV, respectively for the (MnO)4 and (MnO)6 clusters. It 
would be interesting to compare the magnetic structure 
of these clusters with the [Mni20i2(CH3COO)i6(H20)4] 
molecular magnet^. The Mn atoms in the inner shell 
(four Mn'^+) and in the outer sheh (eight Mn^+) are FM 
coupled within both of the shells, but are AFM coupled 
between the shells. This is strikingly different from the 
present (MnO)^^ clusters in the gas phase. 

Similar to the morphological stability, the magnetic 
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stability can also be explained in terms of the total hy- 
bridization index [Fig. [2] (b)] and projected density of 
states (Fig.[3|. Compared to the FM structure, the AFM 
structure has larger orbital spread, and also accompanied 
by higher hybridization index. For example, the 3D-FM 
structure for (MnOje has sharp d-states (Fig. |3|, and 
in contrast, the d-states for 3D- AFM (MnO)6 are more 
widely spread in energy and show larger overlap with 
0-p states. Calculated total hybridization index as 
shown in Fig. [2|b), is much higher in the AFM state 
(8.04) than in the FM state (6.40) in 3D, confirming the 
preference of AFM coupling. We restudied the (MnO)4 
and (MnO)6 clusters using hybrid PBEO functional for 
the exchange-correlation, and we should point out that 
both the structural and magnetic trends discussed here 
are found to be unaltered. 

Compared to AFM bulk MnO crystal^, the overall ex- 
change mechanism is complicated in MnO clusters. The 
semi-empirical Goodenough-Kanamori rules can be em- 
ployed to understand the Mn-Mn magnetic coupling. In 
addition to the Mn-O-Mn superexchange mechanism, di- 
rect Mn-Mn exchange mechanism is also present in these 
clusters since the Mn-O-Mn angle is much smaller than 
180° (Fig. [1]). As we have discussed earlier, pure Mn4 
is FM with AE{FM - AFM) = -0.2 eV, i.e., direct ex- 
change prefers FM coupling. When O atoms are intro- 
duced, the magnetic structure changes to AFM due to 
stronger superexchange coupling that prefers AFM Mn- 
Mn ordering. In contrast, the direct exchange in pure 
Mne is already AFM, and is further stabilized due to the 



AFM superexchange when O is introduced. Similar to 
the ferrimagnetic Fe406 cluster^, these stoichiometric 
MnO clusters also have a very large magnetic exchange 
(at least 0.11 eV per MnO unit), which is much larger 
than the Mn-based molecular magnets^ ^. Therefore, the 
Curie temperature of these (MnO)^^ clusters would be 
much higher than the corresponding AFM bulk MnO (~ 
118 K). This can be exploited to tailor new materials. 

In summary, we have demonstrated, through a rigorous 
and unbiased potential energy search, that the smaller 
stoichiometric MnO clusters show unusual 2D structures, 
and that Mn atoms are AFM coupled. Both these fea- 
tures are explained in terms of the inherent electronic 
structure of these clusters. Present results deviate from 
the earlier theoretical predictions as those works explored 
only a small subspace of the potential energy surface con- 
strained by the high symmetry structures and ferromag- 
netic coupling^^^'^. Although the experimental results on 
such clusters are scarce, the present results agree well 
with the limited experimental predictions on the cluster 
structure and stabilit}^^^. However, there is no experi- 
mental evidence on the evolution of magnetic structure; 
we believe the complementary infrared dissociation spec- 
troscop}/^^ will be helpful to confirm both the geometric 
and magnetic structures of such transition metal oxide 
clusters in the gas phase. 
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